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The halo nucleus 11 Li is treated as a three-body system consisting of an inert core of 9 Li plus 
two valence neutrons. The Faddeev equations are solved using separable potentials to describe 
the two-body interactions, corresponding in the n- 9 Li subsystem to a P1/2 resonance plus a virtual 
s-wave state. The experimental 11 Li energy is taken as input and the 9 Li transverse momentum 
0^ ' distribution in 11 Li is studied. 
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Recent experiments with radioactive beams have permitted the study of properties of nuclei close to the neutron 
drip line, that is, close to the stability line for decay through neutron emission. These nuclei have been shown to 
I/"} ' possess neutron halos characterized by exceptionally large radii and narrow momentum distributions of the decay 
fragments in break-up experiments [1-2]. Among these nuclei special attention has been paid to 11 Li, a nucleus with 
a radius of about 3 fm and a two-neutron separation energy of only 0.3 MeV. In the case of n Li grossly two lines of 
theoretical calculations have been followed; in one, the conventional shell model or Hartree-Fock approach is used [3], 
in the other a cluster model assuming a core of 9 Li plus two neutrons is taken[4-8]. Because of the halo, shell model 
calculations require a very large single particle basis for the diagonalization of the hamiltonian. The cluster model 
seems particularly suited to the case of 11 Li, considering the small 2n separation energy, the fact that 9 Li is a normal 
nucleus with a neutron separation energy of 4 MeV and that 11 Li is a Borromean nucleus, that is, no two-body 
Q\ , subsystem of the three-body system does form a bound state. In the three-body model for 11 Li, calculations are 
0^ " hampered by a lack of information with regard to the n- 9 Li interaction. Several calculations [4- 7, 24] were performed 
, using local potentials for this interaction and there is also one calculation using separable potentials [8]. The parameters 
of the potentials were adjusted to produce a n- 9 Li resonance, which has most frequently been assumed in the P1/2 
i— I | channel, although the experimental data are not conclusive. According to Wilcox et al.[9] a resonance occurs at 
0.80 ± 0.25 MeV. Following more recent work[10] there is a 2 + resonance situated at 0.42 MeV and a 1 + state at 
0.80 MeV. In addition one has now strong evidence for an enhancement of the production of 10 Li near threshold in 
reactions involving 11 Be, n Li and 11 B[11,12]. This is interpreted as due to an intruder virtual s-wave state of the n- 9 Li 
system near zero energy corresponding to a scattering length of -20fm or less [12]. A three-body study [13] indicated 
that this barely unbound state in 10 Li is able to explain the extra narrowness of the momentum distribution of 9 Li in 
the fragmentation of 11 Li. With a s-wave scattering length of -44 fm and a p]y 2 resonance energy of 0.35 MeV, the 
2n separation energy in n Li and the momentum distribution of the 9 Li fragment are fitted in Ref.[13] assuming the 
sudden approximation [14] and neglecting the final state interaction in the break-up reaction. It is our purpose here, 
to make a similar fit applying the three-body model with separable potentials developed earlier[15] and which proved 
to be effective in the description of the structure of 18 and 18 F as well as 16 0(d,p) stripping. 

We consider the halo nucleus 11 Li as a three-body system consisting of a 9 Li core (particle 3), which stays inert, 
plus two valence neutrons (particles 1 and 2). The valence neutrons couple with the orbital angular momentum of 
the core to total angular momentum and parity 3 7T —Q + , the spin and parity of 11 Li being then due to the value | of 
the 9 Li core. The neutron- 9 Li system is assumed to have a pi/ 2 resonance of width T = 0.15±0.07 MeV at an energy 
E r =0.42±0.05 MeV [10] and also a S1/2 virtual state close to zero energy [12,16]. We must also take into account the 
Pauli principle which does not allow a valence neutron to occupy the lsi/2 and the IP3/2 single particle states which 
are already filled in the core. To fix the energies ei Sl/2 and £i P3/2 of these states, we proceed as follows. First by 
doing an interpolation between the experimental value -4.053 MeV of ei P3/2 for A=9 [17] and the values in Fig. 2-30 
of Bohr-Mottelson's book [18] we obtain ei P3/2 ~ —7 MeV for the A=10 system. To obtain ei Sl/2 we suppose that 
the separation between the lsi/ 2 level and the centroid of the levels IP3/2 and lpi/ 2 is Tiuj and use the prescription 
hcu = 45 A~ 1/3 - 25 A~ 2/3 MeV, which is appropriate for light nuclei [19]. For A=10, one has ?icj=15.501 MeV 
and ei Sl/2 results equal to -20.028 MeV. To account for the Pauli blocking of the states ls!/ 2 and IP3/2, we use the 
projection method of Kukulin [20]. From now on, we consider h=l. 
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To describe the neutron- 9 Li interaction we use a separable potential which acts on the S1/2, P3/2 and p x / 2 waves: 

lj a 

* £ <Pii Wj >> <w,> ip;> (* = l > 2 ) . (!) 

where lj = 0|, 1 § and 1 i , m is the reduced mass of the n- 9 Li system (m=^M, M being the nucleon mass) and 
Pi is the momentum of neutron i with respect to 9 Li. The Si/ 2 potential is a three-term (a=l,2,3) potential, with 
the form factors chosen as 

= (« 2 + "oi) exp(-/? V/2) , (2) 
«g(<?) = (^-/3 V)exp(-/3 V/2) , (3) 

v^(q) = exp(-/3 V/2) . (4) 

The first term is chosen in such a way that, alone, it reproduces the lsi/2 state. Taking for a i and A^V the special 
values 



a i = J2m\e u | (5) 
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we get a bound state of energy ei Sl , 2 and wave function 

*o^(q) = N oi exp (-/3 V/2) tf i M (q) , (7) 

which is precisely the lsi/ 2 oscillator function in momentum space if we set /3o = l/y/mu) . From the values chosen 

for ei Sl/2 and w, we obtain a i=0.932 fm" 1 , A^V =8.417 fm 5 and /3 = 1.724 fm. 

The addition of the second term to the Si/ 2 potential does not affect the bound state generated by the first term. 
This is a consequence of the orthogonality relation 
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dqq 2 [exp (-/3 V/2)] [Q - /? V) exp (-(3 2 q 2 /2)} = . (8) 



However, the scattering states are affected and the two terms together can give rise to a virtual S1/2 state. For 
(2) 

A^7 =2.696 fm, we obtain a virtual state placed at an energy e v —-40 keV on the second Riemann sheet. The cor- 
u 2 

responding scattering length is a Sl/2 ~-20 fm. We remark here that the virtual state has a 2si/ 2 character, since by 
increasing slightly the strength of the second term (a=2) in the Si/ 2 wave potential, it becomes a bound state with 
two nodes. 

The third term of the Si/ 2 potential is the projection operator for the forbidden lS]y 2 state, constructed according 
to the prescription of Kukulin[20]. The corresponding form factor (Eq.(4)), being proportional to the wave function 
of the bound state produced by the first two terms, is orthogonal to the scattering states generated by these terms. 
Therefore the third term does not affect the scattering and the virtual state remains unchanged. It is not so for the 
lsi/ 2 bound state. Although the wave function given by expression (7) remains an eigenfunction for the three-term 

potential, the corresponding energy is affected. By considering the third term repulsive (A^ < ), we remove the 

u 2 

bound state to the continuum part of the spectrum (therefore, it becomes a continuum bound state) and, by making 



2 



— > — oo, the forbidden lsi/ 2 state is projected out. 
For the p 3 / 2 potential, we consider a two-term potential with form factors 

(9) = 5 (<Z 2 + ^ 3 ) cxp /2) , (9) 

^(g) = qexp(-0*<?/2) . (10) 



The first term of the potential is chosen so as to reproduce the IP3/2 bound state. With the choice a 1 3 = 2m |ei P3/2 1 
and 
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the potential produces a bound state of energy e\ P3/2 and wave function identical to the lp 3 / 2 harmonic oscillator 
wave function 

= 7V 1§ q exp (-0tf/2) ^(q) . (12) 
For the choice /?i=/3 =1.724 fm and ei p =-7 MeV, we obtain a l2 =0.551 fin" 1 and A^V =20.018 fm 7 . 

(21 

The second term in the p 3 / 2 potential, with A^ — > —00, is nothing but the projection operator which projects out 

the forbidden lp 3 / 2 state, the scattering states remaining the same as those produced by the first term alone. Thus, 
the p 3 / 2 phase shift is dominated by the occupied lp 3 / 2 bound state, there being no resonances as it might occur in 
the case of a local potential. 

The p!/2 potential is taken as a one term potential with the form factor given by 

(q) =q{q 2 + a\ h ) cxp (-/3 2 g 2 /2) • (13) 

Thus, t> is of the same form as v^}. For simplicity we have also assumed the same parameter /3i in the exponential 
1 2 1 2 

part of the pi/ 2 and p 3 / 2 form factors. Wc initially make the choice a x i = a^a, thus getting 0^1=0.551 fm -1 . Using 

this value and the condition that the p!/ 2 resonance occurs at 0.42 MeV, we obtain =13.535 fm 7 . The width of 

1 2 

the resonance turns out to be 0.12 MeV and compares with the experimental value mentioned before. 

Regarding the interaction between the valence neutrons, we assume a free nucleon-nucleon interaction considering 
the low nucleon density in the region of the halo. For the + state which we are considering, if one assumes pure single 
particle harmonic oscillator states, the (pi/ 2 ) 2 configuration is a superposition of both the 1 So and 3 Pi states (33 and 
67 percent respectively) of the n-n subsystem, while only -"-So appears in the (s!/ 2 ) 2 configuration. Calculations by 
other authors (for instance, Ref.[13]) have shown that the 3 Pi potential, which is repulsive, changes the 2n separation 
energy by at least 50%. We therefore use a potential which acts in both 1 So and 3 Pi channels, and take it separable, 

(p|Vi 2 |p') = J2 ^"^p VXSI W v *si(p') 

XSI 

*E wvls 1 ) Ms 1 m . a 4 ) 

Mi 

where p is the relative momentum between the two neutrons, (A5'/)=(000) and (111) and the form factors are of the 
Yamaguchi type: 

«ooo(p) = [p 2 + aooo] , (15) 

viu{p) = V [p 2 + "111] 2 • (16) 

Using the values ais =-17 fm and rig =2.84 fm for the neutron-neutron scattering length and effective range in the 
1 Sq channel [21] and a3 Pl =2.2 fm 3 and r3 Pl =-8.0 fm -1 for the 3 Pi channel [22], we fix the parameters of the V12 
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potential as A oo = 1-662 fm -3 , a oo — 1-130 fm -1 , A m = —0.078 fm~ 5 and am = 0.693 fm -1 . The negative 
value of Am means that the 3 Pi term is repulsive. 

The proposed interactions are used as input in the homogeneous Faddeev equations 

= Go Ti (V 2) + * (3) ) , 
* (2) = Go T 2 (V 3 > + , (17) 

* (3) = Go lb + * (2) ) , 

where , and ~$>^ are the Faddeev components of the total wave function, 

$ = + *(2) + *(3) , (18) 

Go = (E — Ho) 1 is the Green's function and T\,T 2 and Tyi are the T-matrices corresponding to the potentials 
V\, V2 and V12 respectively. After performing the angular momentum decomposition, we end up with a homogeneous 
system of coupled integral equations in one variable. The equations are then transformed into a set of homogeneous 
algebraic equations using the Gauss quadrature method to approximate each integral by a finite sum. The zero of 
the determinant of this system of equations gives the separation energy, S 2n , of the two valence neutrons. From 
the corresponding three-body wave function, we calculate the momentum distributions. We restrict ourselves to 
the calculation of the transverse momentum distribution of the 9 Li core in n Li. We should mention that, in our 
model, the transverse and the parallel momentum distributions turn out to be identical. An approximate equality has 
indeed been verified experimentally for the 9 Li fragment in the break-up of n Li [23]. Corrections due to final state 
interactions in the case of the momentum distribution of the core fragment should be small according to reference 
[24]. 

For the chosen values of the parameters, we obtain the value S 2n =0.293 MeV which agrees with the (average) 
experimental value 0.294±0.030 MeV, reported in reference [25]. However, the width of the transverse momentum 
distribution of 9 Li turns out to be too large, as shown in the figure (dotted line). The experimental data are from 
Kobayashi et al. [26] (actually, the experimental points shown in the figure are taken from Fig. 2 of reference [6]). We 
here mention that, if the 3 Pi component of the n-n interaction is suppressed the value of S 2n duplicates, becoming 
equal to 0.597 MeV. 

The result may be improved by taking for (3q a value smaller than the one given by the prescription (3$ — 1/ ' \fmu) 
= 1.724fm. We made a new calculation in which the S1/2 and P1/2 potentials are modified as follows. Assuming 

/3o=l-5fm and maintaining the previous values ei Sl/2 = -20.028MeV and a o i=0.932 fm -1 , we determine to be 

/ 2 U, 2 

4.959 fm 5 . By setting = 2.219 fm the virtual state is positioned at e„=-40keV, the corresponding scattering 

length being a Sl/2 =-20fm. To determine the pi/ 2 potential, we keep (3\ unchanged and adjust and a 1 1 with the 

condition that the calculated p 1 / 2 resonance energy occurs at 0.42 MeV, the resonance width T remains inside the 
range 0.15±0.07 MeV and, in addition, that the calculated 2n separation energy agrees with the experimental value 
0.294MeV. The parameters of the p 3 / 2 interaction are the same as in the previous calculation. The corresponding 
transverse momentum distribution is given by the dashed line of the figure. The half width at half maximum(HWHM) 
is 40 MeV/c and the three-body state is a superposition mainly of 1 So (47%) and 3 Pi (48%) configurations of the n-n 
subsystem. 

Repeating the previous calculation, with a still smaller value, /3o=l-4 fm, we get 35 MeV/c for the HWHM. The 
agreement with experiment is very good (solid line of the figure). A similar result for the transverse momentum 
distribution is also obtained taking (3q=1.5 fm, and changing the scattering length to a Sl/2 =-40 fm (dot-dashed line). 
In Table 1 we summarize the parameters of the two-body Si/ 2 and p!/ 2 potentials used in the curves I-IV and in 
Table 2 the corresponding three-body results. The increase of the n-n 1 S contribution in the cases II-IV, points to 
an increase of the contribution of the (si/2) 2 configuration to the three-body wave function. 

Our calculation shows that relevant 11 Li data may be fitted by a three-body model using simple separable potentials 
with parameters adjusted to two body data. The fact that one needs for the parameter (3 a value not in accordance 
with the prescription (3q = 1/ ' ^Jmui indicates that the effective neutron- 9 Li potential may differ appreciably from 
the usual single-particle potential. The use of a local neutron- 9 Li potential should corroborate this conclusion. In 
fact, one can verify, by using for instance a simple square-well potential, that in order to obtain an intruder 2si/2 
state near zero energy a much deeper potential than usually is needed. Such odd potentials may indicate that a more 
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detailed description, which does not consider the core as a structureless object, is required. 
The numerical calculations were performed at LCCA-USP. 
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TABLE I. Parameters of the two-body s 1( / 2 and p 1(/2 potential corresponding to the curves I- IV presented in Fig. 1. These 
parameters correspond to eis 1/2 = -20.028MeV and E r =0.42MeV. For the p 3 / 2 potential we use /?i=1.724fm, A (1 ^=20.018fm" 1 , 

/ 1, 2 

a 1 3=0.551fm _1 , corresponding to ei P3/2 =-7MeV. The width T of the p 1 / 2 resonance and the scattering length a Sl/2 are also 



given. 
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TABLE II. Three-body results corresponding to the parameters listed in Table 1. S 2 „ is the 2n separation energy from 11 Li 
and 7 is the HWHM of the calculated momentum distribution of the 9 Li core in 11 Li. The last four columns give the fractional 
admixture of n-n states in the lx Li wave function. 
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FIG. 1. Transverse momentum distribution of 9 Li in the 11 Li. The squares and circles are experimental data[26] corresponding 
to Px<0 and P^>0 respectively. Curve (I) corresponds to /3o=/3i=1.724fm and scattering length a Sl/2 =-20fm, (II) to /3o = 1.5fm, 
/ 9i=1.724fm, a Sl/2 =-20fm, (III) to /3 =1.4fm, /3i=1.724fm, a Sl/2 =-20fm and (IV) to /3 =1.5fm, /3i=1.724fm, a Sl/2 =-40fm. 
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